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ABSTRACT 



We describe a method by which gravitational wave observations of eccentric 
binary systems could be used to test General Relativity's prediction that gravi- 
tational waves are dispersionless. We present our results in terms of the graviton 
having a non-zero rest mass, or equivalently a non-infinite Compton wavelength. 
We make a rough estimate of the bounds that might be obtained following grav- 
itational wave detections by the space-based LISA interferometer. The bounds 
we find are comparable to those obtainable from a method proposed by Will, 
and several orders of magnitude stronger than other dynamic (i.e. gravitational 
wave based) tests that have been proposed. The method described here has the 
advantage over those proposed previously of being simple to apply, as it does not 
require the inspiral to be in the strong field regime nor correlation with electro- 
magnetic signals. We compare our results with those obtained from static (i.e. 
non-gravitational wave based) tests. 

Subject headings: gravitation — gravitational waves — relativity 



In Einstein's theory of General Relativity, linearization of the field equations shows that 
small perturbations of the metric obey a wave equation (Misner, Thorne & Wheeler 1973). 
These small disturbances, referred to as gravitational waves, travel at the speed of light. 
However, some other gravity theories predict a dispersive propagation (see Will & Yunes 
(2004) for references). The most commonly considered form of dispersion supposes that the 
waves obey a Klein-Gordan type equation: 
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Physically, the dispersive term is ascribed to the quantum of gravitation having a non-zero 
rest mass m g , or equivalently a non-infinite Compton wavelength A g = h/m g c. The group 
velocity of propagation for a wave of frequency / gw is then 



valid for A gw <C A s ; only in the infinite frequency limit is General Relativity recovered, with 
waves traveling at the speed of light (Will 1998). 

Over the past few decades a number of different dynamic tests of this dispersive hy- 
pothesis have been described, i.e. tests making use of direct observations of gravitational 
waves or their radiation reaction effects (Cutler et al. 2003; Finn & Sutton 2002; Larson & 
Hiscock 2000; Will 1998; Will & Yunes 2004). 

In this paper we add another method to this list; we consider gravitational radiation 
from eccentric binary systems. Such binaries emit gravitational radiation at (infinitely many) 
harmonics of the orbital frequency (Peters & Mathews 1963). Our idea lies simply in mea- 
suring the phase of arrival of these harmonics. Dispersion of the form described by equation 
(2) would be signaled by the higher harmonics arriving slightly earlier than the lower har- 
monics, as compared to the General Relativistic waveform. We present a rough estimate of 
the bounds that might be obtained, deferring a more accurate calculation to a future study 
(Barack & Jones, in preparation). 

The plan of this paper is as follows. In §2 we derive formulae to make a simple estimate 
of the bounds that might be obtained using our method. In §3 we estimate bounds obtainable 
on A g for LISA observations of two sorts of binary systems. Finally in §4 we summarize our 
findings and compare with those of other authors. 



To derive a rigorous estimate of the bound one should add the graviton mass to the 
list of unknown source parameters to be extracted from the measured signal, as was done 
by Will in the case of circular orbits (Will 1998). The m g -dependent waveform can then 
be computed, allowing calculation of the Fisher information matrix T a i,, which could then 
be inverted, the r~* mg component, evaluated at m g = 0, giving the best bound obtainable 
(Will 1998). For the case of eccentric binaries such a calculation is not easy, and so in this 
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Derivation of the bound 
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General formula 
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paper we make a preliminary estimate of the possible bounds, without going to the trouble 
of calculating T ab . 

We will begin by deriving a general formula for estimating the bound on A g that could 
be obtained from a system which produces gravitational waves at two different frequencies, 
say /gw,i and /gw,2- The two gravitational waves will travel with (different) speeds v gt i and 
v gi 2, and so their journey times to the detector a distance d away will differ by a time interval 
At given by 

^ r 1 11 

(3) 
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'gw,l Jgw,2. 

Multiplying this by 2nf gw , where f gw is a characteristic frequency in the problem, gives 
the accumulated difference in phase of arrival of the two signals caused by the dispersion, 
measured in terms of radians of phase of f gw : 
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This is to be compared with the accuracy with which the phase of arrival of the waves can 
be extracted from the noisy gravitational wave data stream. In the high signal to noise ratio 
regime the error in extracting the phase of a continuous signal can be written as 
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where we follow the notation of Cutler et al. (2003). In this formula p is the signal to noise 
ratio of the measurement and a is a dimensionless factor that depends upon how many 
unknown parameters (including the phase) need be extracted from the signal. 



The lower bound that can be placed on A g comes from equating A$ dispcrsion and A$ c 
to give: 

pd 



a 



P 



P 



(6) 



.Jgw,l Jgw,2. 

This shows that the best bounds will come from high mass (i.e. high pd), high eccentricity, 
low orbital frequency systems. 



2.2. Application to eccentric binary systems 

We will now apply this method of estimation to eccentric binary systems. In general 
many more than two harmonics will contribute significantly to p, so equation (6) is not 
directly applicable. In order to take advantage of this spread we will make the following 
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identifications. We will set p equal to the total signal to noise of the observation. To identify 
appropriate frequencies, consider a plot of the signal to noise of the n-th harmonic, p n , verses 
gravitational wave frequency n/ or bit- We will set f gw to be the frequency at which this curve 
peaks, and / gw ,i,/gw,2 as the frequencies corresponding to the lower and upper full-width- 
at-half-maximum. In reality only discrete harmonic frequencies exist, but for the purpose 
of defining f\, f 2 and / gw , we will treat the curve as continuous, interpolating to find the 
necessary frequencies. (A formalism using only discrete frequencies would have introduced 
spurious step-wise changes in our bounds on A g as a function of eccentricity). 

Identification of a suitable a value, which quantifies the error in phase measurement, 
is more problematic. Cutler et al. (2003) consider errors in measuring the phase of a single 
monochromatic signal of known sky location; they find that a < 3 for a large fraction of the 
possible binary orientations. Barack & Cutler (2004) examine extreme mass ratio inspirals. 
They find phase measurement errors which again yield a ~ 3 (see the A(to)v parameter of 
their Table III). 

However, even in the dispersionless case of General Relativity, the relative phasing of 
the detected harmonics is non-trivially determined by the source's sky location and the 
relative orientation of the detector and binary system (Barack & Cutler 2004). The phase 
differences we are considering here are additional delays caused by dispersive propagation. 
Clearly, then, the results of Cutler et al. (2003) and Barack & Cutler (2004) do not directly 
apply to our problem. Only a full Fisher matrix calculation will accurately show how we can 
disentangle the phase differences contributed by dispersion, measurement error and those 
intrinsic to the binary. We expect that in those situations where the system parameters, 
including its sky location and orientation relative to LISA, are measured accurately, the 
dispersion-induced phase delays will be measured accurately too. In the absence of a full 
Fisher matrix calculation to evaluate the correct measurement errors we will set a = 10, but 
note that this is the weakest link in our estimate. 

In particular, if the various geometric factors that enter the problem conspire such that 
a dispersionless signal from a certain binary is very similar to the dispersed signal from a 
binary with slightly different parameters (e.g. a slightly different sky location), then the 
errors in A$ could be very much larger than estimated here. Also, a will depend upon the 
type of system being studied. It will generally be smaller for systems where information in 
addition to the gravitational wave signal is available, e.g. Galactic binary systems where 
optical measurements give accurate sky locations. Note, however, that a enters the bound 
on A g only rather weakly, as A g oc and so we hope that our ignorance of this factor 

will not change our qualitative conclusions. 
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3. Results 

It is expected that gravitational radiation reaction will result in most binary systems 
detectable by ground based interferometers being nearly perfectly circular (Peters 1964) 
and so will be unusable for deriving a bound of the sort described here. We will therefore 
concentrate exclusively on (two sorts of) binaries in the LISA band. In equation (6) we will 
set a = 10, as discussed above. When calculating p we will assume an integration time of one 
year. We computed the LISA noise using the Online Sensitivity Curve Generator 1 (Larson 
2003), which included a fit to the Galactic white dwarf background (Bender & Hils 1997). 

3.1. Extra-Galactic extreme mass-ratio binaries 

We consider here the inspiral of a solar-mass type black hole into a massive one. These 
are excellent systems from our point of view, as they are expected to dominate the LISA 
inspiral event rate and, crucially, many will have very large eccentricities (Barack & Cutler 
2004; Gair et al. 2004). 

To see if such systems can indeed be used to obtain a bound on A g , in Figure 1 we plot the 
eccentricity-orbital frequency phase space for a (10 6 , 1O 1 )M binary at a distance of 1 Gpc. 
The upper curve describes the innermost stable orbit (ISO) (Will 1998); binary systems in 
Nature only exist below this curve. The lower curve gives the minimum eccentricity required 
for the system to be detectable, with multiple harmonics contributing significantly to p. 
[Our exact criterion is to see if p exceeds some detection threshold p m j n when the single 
strongest harmonic is removed from the sum. We have set p m i n = 15, as would be reasonable 
if computational power does not limit the search (Gair et al. 2004)]. Our methods are only 
applicable for systems above this curve. It follows that we can use binary systems which lie 
between these two curves to bound A g . Fortunately we see that this means that binaries in 
a significant portion of the e — / plane are of use to us. To illustrate this, a trajectory of a 
plausible LISA source is shown between the two curves, with an eccentricity at the ISO of 
about 0.24. This system spends about 10 years between the two curves. 

In Figure 2 we show the actual bounds on A g that could be obtained from observations 
of extreme mass ratio systems. The distance is still fixed at 1 Gpc, but now we fix the 
orbital frequency at 10~ 3 Hz and leave the eccentricity as a free parameter. Results for 
binary systems with Mi = 1O 6 M and several different values of M 2 are given, as indicated. 



http://www.srl.caltech.edu/~shane/sensitivity 
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The following features are of note: (i) Each curve terminates at a minimum eccentricity 
below which the system is undetectable and/or fewer than two harmonics contribute signifi- 
cantly to p, and at a maximum eccentricity above which the system is dynamically unstable, 
(ii) For a system of given masses, the bound increases slightly (i.e. becomes stronger) the 
larger the eccentricity, (iii) Stronger bounds are obtained from more massive systems, and 
can be obtained for wider ranges of the eccentricity. 

3.2. Stellar mass Galactic binary systems 

LISA will be able to detect gravitational waves from a large number of low mass Galactic 
binaries, consisting of white dwarfs and/or neutron stars (Danzmann et al. 1996). To inves- 
tigate the suitability of these systems for bounding A g , in Figure 3 we plot the eccentricity- 
frequency phase space for a Galactic (1, 1)M Q binary at a distance lkpc. We set p m i n = 8, 
although a lower value could be used for electromagnetically studied binaries (Danzmann 
et al. 1996). We do not show the ISO curve as for all plausible eccentricities such a binary 
would go dynamically unstable in the much higher LIGO frequency band. Clearly, binaries 
in a large portion of the phase space are of use for bounding A g . However, unlike the case of 
the extreme mass ratio inspiral, there is no compelling reason to expect the eccentricities of 
these systems to be large. Many of them will have gone through a period of mass transfer in 
the past, which is believed to be an efficient circularizer. Nevertheless, as we require merely 
one or more of them to have a sufficiently large eccentricity, greater than about 5 x 10~ 3 , a 
bound on A g may well be obtained. 

In Figure 4 we present the bounds on A g that would be obtained from observations of 
various equal-mass binaries at a distance of 1 kpc and with an orbital frequency 10~ 3 Hz. 
The qualitative form is the same as in Figure 2, except we terminate the curves at the high 
eccentricity end at e = 0.55 as such extreme eccentricities seem unlikely. 

4. Comparison with previous methods and summary 

In Table 1 we collect together reported and proposed dynamic bounds on A g that have 
appeared in the literature, and add two proposed bounds from this work. As is clear from 
perusal of the Table and Figures 2 and 4, the bounds presented here for low mass Galactic 
systems are comparable to those of Cutler et al. (2003), while our bounds from extreme 
mass ratio inspirals are comparable to those of Will & Yunes (2004) for massive black hole 
coalescence. It should be remembered that our numbers can only be regarded as estimates, 
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particularly given our rough guess as to the accuracy with which dispersion-induced phase 
delays can be measured. However, even if a, the parameter which quantifies this error, 
were four orders of magnitude larger than assumed here, our results for extreme mass ratio 
inspirals would still beat both solar system and Galactic low mass binary tests. 

The method presented here has several advantages over other methods. The analysis 
of Cutler et al. (2003) requires knowledge of the initial relative phases of the X-ray and 
gravitational wave signals from an accreting white dwarf system; it is not clear if the accretion 
process will be sufficiently well understood to allow this. Less problematically, the method of 
Will (1998) requires knowledge of the phasing of the binary inspiral waveform in the strongly 
chirping regime, as it is this frequency variation that allows the dispersion test. In contrast, 
the method described here is very simple, requiring only that multiple harmonics can be 
detected. It is not necessary for the binary to be chirping significantly, and correlation with 
other (i.e. non-gravitational) radiation is not required. 

Returning to equation (1), in the static regime the solution is of the form of a Yukawa- 
type potential, i.e. a Newtonian one suppressed by an exponential exp (— r/A g ). This offers 
the possibility of bounding A g by looking for departures from Newtonian gravity in the non- 
radiative regime. Such results are given in Table 2. Talmadge et al. (1988) used planetary 
ephemeris data to obtain their bound, while Goldhaber & Nieto (1974) cited evidence of 
gravitational binding of galaxy clusters, suggesting that the exponential suppression is not 
important over length scales of the order of a Mpc. The bounds that could be obtained by 
using the methods described in this paper would be better than the solar system bounds by 
around 5 orders of magnitude. However, they are weaker than those from galaxy clusters by 
3 orders of magnitude. 

Therefore, if equation (1) is the correct linearization of the true theory of gravity, and if 
galaxy clusters are indeed gravitationally bound, the bounds on X g from the dynamic sector 
are much weaker than those from the static sector. However, the possibility remains that 
equation (1) is not the correct linearization, the static potential is not suppressed, but the 
wave propagation is nevertheless dispersive, i.e. equation (2) holds but is not derived from 
an equation of the form of equation (1). The only way of settling this is to use the methods 
proposed in the dynamic regime. It could even be the case that neither equations (1) nor (2) 
are correct, but that gravitational waves have some other form of dispersion. The method 
considered here (or any of the methods referred to in Table 1) could be used to identify this. 

To sum up, the estimates in this paper indicate that our method will give bounds on A g 
which are stronger than most other dynamic bounds, being rivaled only by those of Will & 
Yunes (2004). 
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Having used simple estimates to establish the competitiveness of the method presented 
here with other dynamic tests, we are currently working to improve the accuracy of our 
calculation by using the Fisher information matrix to calculate the bound (rather than the 
methods of §2; Barack & Jones, in preparation). We also aim to extend the scope of the 
investigation by considering the full range of anticipated gravitational wave sources for both 
ground and space-based detectors. 

It is a pleasure to thank Leor Barack, Shane Larson, Ben Owen, Steinn Sigurdsson and 
Nico Yunes for useful discussions during this investigation, and the anonymous referee for 
providing comments which improved the manuscript. The Center for Gravitational Wave 
Physics is supported by the National Science Foundation under cooperative agreement PHY 
01-14375. 
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Table 1. Actual/proposed bounds on A g from the dynamic sector 



Reference 


Binary system 


A g bound 
(m) 


1 


Radio pulsars 


1.6 x 10 13 


2 


4U1820-30 


1 x 10 16 


3 


e = 0.1, (O.5,O.5)M 


1 x 10 16 


2 


Ideal low mass binary 


1 x 10 17 


4 


(1O 4 -1O 4 )M 


4 x 10 18 


3 


e = 0.3, (1O 6 ,1O 2 )M 


3 x 10 19 


4 


(1O 7 -1O 7 )M 


5 x 10 19 



References. - (1) Finn & Sutton (2002); (2) 
Cutler et al. (2003); (3) This work; (4) Will & 
Yunes (2004). 
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Table 2. Actual bounds on Ag from the static sector 



Reference 


System 


A g bound 






(m) 


1 


Solar system 


2.8 x 10 15 


2 


Galaxy clusters 


~ 10 23 



References. (1) Talmadge et al. 

(1988); (2) Goldhaber & Nieto (1974). 
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Fig. 1. — The eccentricity-orbital frequency phase space for a (10 6 , 1O)M binary at a 
distance of 1 Gpc. The ISO curve, a sample trajectory, and the minimum eccentricity (p m i n = 
15) curve are shown. 

Fig. 2. — Bounds on A g obtainable from extreme mass ratio binaries with /orbit = 10~ 3 Hz, 
Mi = 1O 6 M , p min = 15 at a distance of 1 Gpc. 

Fig. 3. — The eccentricity-orbital frequency phase space for a (1, 1)M Q binary at a distance 
of 1 kpc. A sample trajectory and the minimum eccentricity curve (p m i n = 8) are shown. 

Fig. 4. — Bounds on A g obtainable from equal low mass Galactic binaries with /orbit = 
lCT 3 Hz, p min = 8 at a distance of lkpc. 
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